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I. INTRODUCTION 

The electric dipole polarizabilities, a, of ions are im- 
portant to determine the optical properties of ionic crys- 
tals. In addition, for closed-shell ions like the singly ion- 
ized alkali atoms, a is a measure of the core-polarization 
effects in the neutral species. It is, however, nontriv- 
ial to measure a of ions. For the singly charged alkali 
ions, an indirect method to determine a is through the 
measurement of the transition energy between the non- 
penetrating Rydberg states of the neutral species [T] and 
it has been used to determine the a of Cs"*" [3 [S] . In ab- 
sence of experimental data, there is a need for accurate 
theoretical calculations. In the case of neutral atoms, 
accurate values of polarizabilities are essential in studies 
related to the parity non-conservation in atoms ^ , opti- 
cal atomic clocks [3 [6] and physics with the condensates 
of dilute atomic gases [7HS] are of current interest. 

Theoretically, methods based on a wide range of atomic 
many-body theories have been used to calculate a. In 
this regard, the recent review [TD] provides description 
about the various theoretical methods used to calculate 
a. In the present work we use the perturbed relativis- 
tic coupled-cluster (PRCC) theory, which was earlier ap- 
plied the noble gas atoms [TTl [T^], to compute the a 
of the singly charged alkali ions. The PRCC theory is 
an extension of the standard relativistic coupled-cluster 
(RCC) theory to include an additional perturbation and 
for this, we introduce a new set of cluster operators. The 
formulation is, however, general enough to incorporate 
any perturbation Hamiltonian. It must be emphasized 
that, compared to other many-body methods, the use of 
PRCC is an attractive option as it is based on coupled- 
cluster theory (CCT) [13l [Mj: an all order many-body 
theory considered to be reliable and powerful. The re- 
cent review |15| provides an overview of CCT, and vari- 
ants of CCT developed for structure and properties cal- 
culations. The theory has been widely used for atomic 
[TBHTO] . molecular [5D], nuclear [2T| and condensed mat- 
ter physics fF3 calculations. Coming back to the PRCC 
theory, it is different from the other RCC based theo- 
ries in a number of ways, but the most important one is 
the representation of the cluster operators. In the PRCC 



theory, the cluster operators can be scalar or rank one 
tensor operators and it is decided based on the nature of 
the perturbation in the electronic sector. Consequently, 
the theory is suitable to incorporate multiple perturba- 
tions of different ranks in the electronic sector. 

One basic advantage of PRCC theory is, it does away 
with the summation over intermediate states in the first 
order time-independent perturbation theory. The sum- 
mation is subsumed in the perturbed cluster amplitudes 
and this offers significant advantages in computing prop- 
erties like a which involves summation over a complete 
set of intermediate states. 

The paper is organized as follows. In the Sec. II, for 
completeness and easy reference we briefly describe the 
RCC and PRCC theories with the Breit interaction. In 
Sec. IV we introduce the formal expression of the dipole 
polarizability and its representation in the PRCC theory. 
In the subsequent sections we describe the calculational 
part, and present the results and discussions. We then 
end with conclusions. All the results presented in this 
work and related calculations are in atomic units ( h = 
TOg = e = 47reo = 1). In this system of units the velocity 
of light is a~^, the inverse of flne structure constant. For 
which we use the value of = 137.035 999 074 [53]. 



II. OVERVIEW OF THE COUPLED-CLUSTER 
THEORY 

The detailed description of the RCC and PRCC theo- 
ries are given in our previous works. However, for com- 
pleteness and easy reference we provide a brief overview 
in this section. 



A. RCC theory 

The Dirac-Coulomb-Breit Hamiltonian, denoted by 
jjBCB ^ jg g^j^ appropriate choice to include the relativistic 
effects in the structure and property calculations of high- 
Z atoms and ions. There are, however, complications 
associated with the negative energy continuum states of 
^DCB xhese lead to variational collapse and continuum 
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dissolution [23]. One remedy to avoid these complica- 
tions is to use the no-virtual pair approximation. In this 
approximation, for a singly charged ion of N electrons 
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where a and /3 are the Dirac matrices, A-|-+ is an op- 
erator which projects to the positive energy solutions 
and VN{ri) is the electrostatic potential arising from 
the Z = {N + 1) nucleus. Projecting the Hamiltonian 
with ensures that the effects of the negative en- 

ergy continuum states are removed from the calculations. 
The last two terms in l/r^j and g^^r^j), are the 

Coulomb and Breit interactions, respectively. The later, 
Breit interaction, represents the transverse photon inter- 
action and is given by 
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The Hamiltonian satisfies the eigen-value equation 
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where, is the exact atomic state. In CCT the exact 
atomic state is defined as 



(4) 



where is the reference state wave-function and T*^*^^ is 
the unperturbed cluster operator, which incorporates the 
residual Coulomb interaction to all orders. We have in- 
troduced the superscript to distinguish it from the second 
set of cluster operators, the perturbed cluster operators, 
to be introduced later. In the case of a closed-shell ion, 
the model space of the ground state consists of a single 



Slater determinant, |$o), and T(°) = E,=i where ^ 
is the order of excitation. However, in actual computa- 
tions, incorporating T^^^ with j ^ 4 is difficult with the 
existing computational resources. A simplified, but quite 
accurate approximation is the coupled-cluster single and 
double (CCSD) excitation approximation, in which 
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This is an approximation which embodies all the impor- 
tant electron correlation effects, and is a good starting 
point for structure and properties calculations of closed- 
shell ions. In the second quantized notations 



(6a) 
(6b) 



where t■^^ are cluster amplitudes, a] (a^) are single parti- 
cle creation (annihilation) operators and abc. . . {pgr . . .) 
represent core (virtual) states. For the present work, the 
ground state is the required atomic state |^'o) = e 
and satisfies the eigenvalue equation 



(7) 



where, Eq and |$o) are the energy and reference state of 
the ground state, respectively. Following similar proce- 
dure, the CC eigenvalue equation of the one-valence ^2S] 
and two-valence [27j systems may be defined. 



B. PRCC Theory 

In the PRCC theory we introduce a new set of cluster 
operators, T^-*-', to incorporate an interaction Hamilto- 
nian, -ffint, perturbatively. For general representation, 
we consider T^-*") as tensor operators of arbitrary rank 
and depends on the multipole structure of Hmt- The 
new cluster operators follow the selection rules associ- 
ated with Hint and the modified ground state eigenvalue 
equation, after including the perturbation, is 



DCB 



Ai/i„t)|*o) =i^o|*o) 
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where A is the perturbation parameter, l^to) is the per- 
turbed ground state and Eq is the corresponding eigen 
energy. To calculate the electric dipole polarizability, a, 
consider the perturbation as the interaction with an elec- 
trostatic field E. The interaction Hamiltonian is then 
ffint = — J2i • E = D • E, where D is the many elec- 
tron electric dipole operator. Here, Hint is odd in parity 
and to be more precise, D, the operator in the electronic 
space is odd in parity and a rank one operator. Hence, 
the cluster opertors T'--'--' are also rank one tensor oper- 
ators and odd in parity, meaning, they connect states of 
different parities. Further more, the first energy correc- 
tion (4'o|^^int|*o) = and therefore, Eq = Eq. We can 
then write, using PRCC theory, the perturbed ground 
state as 



l*o) 



l^o) 



1 + AT(^' • E 



where, we have introduced the scalar product between 
T^^^ and E for a consistent representation of the states 
and operators. The advatage of introducing T^-^^ and 
using l^'o) is, it allows a systematic consolidation of the 
correlation effects arising from multiple perturbations. 

Based on the analysis of the low-order many-body per- 
turbation theory diagrams, the single and double excita- 
tion operators of PRCC theory are represented as 



(10a) 
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where T^^^ are the cluster amphtudes and Ci(r) are C- 
tensors of rank i. To represent T^^\ a rank one operator, 
we have used the C-tensor of similar rank Ci(f). And, 
the key difference of T^^-* from T^^'^ is la + lp must be odd, 
in other words (— l)'»+'p — —1, where, la {Ip) is the or- 
bital angular momentum of the core (virtual) state a (p) . 
Coming to T2^'', to represent it two C-tcnsor operators of 
rank I and k are coupled to a rank one tensor operator. In 
terms of selection rules, the angular momenta of the or- 
bitals and multipoles in T2^'' must satisfy the triangular 

conditions \ ja-jp\ < I ^ (ja+jp), \jb-jq\ < < (jb+jq) 

and |Z — fc| ^ 1 ^ {I + k). The other selection rule fol- 
lows from the parity of iJint , the orbital angular momenta 
must satisfy the condition (— l)('»+'p) = _(_i)('i.+'<j). 



C. PRCC equations 

The ground state eigenvalue equation, in terms of the 
PRCC state, is 

In the CCSD approximation we define the perturbed 
cluster operator T^^^ as 



^ rj,(l) 



.(1) 
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Using this, the PRCC equations are derived from Eq. 
(111. The derivation involves several operator contrac- 



tions and these are more transparent with the normal 
ordered Hamiltonian H^'^^ = i/^CB _ ($^|£fDCB|^.^ 
The eigenvalue equation then assumes the form 



[i/DCB + Ai?i„t] |*o> - [Eo - ($o|i^°^^l<i>o)] l*o)- 

(13) 

A more convenient form of the equation is 
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where, AEo ^ Eq ~ {^o\H^'^^\^q) is the ground state 
correlation energy. Following the definition in Eq. (|9]), 
the PRCC eigen-value equation is 

(15) 



Applying e 



_rp{0) 



from the left, we get 



where i/DCB ^ 



.AT 



_rp(0-) 



'^'"'a>o)=Ai?oe^T'^'-E|$o), (16) 
^DCBgT*"' jg ^j^g similarity trans- 



formed Hamiltonian. Multiplying Eq. ( 16 1 from left by 
e"^"^' and considering the terms linear in A, we get the 
PRCC equation 



fVDC rp(l) 



E-Hffi„t|$o) =0. 
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Here, the similarity transformed interaction Hamiltonian 
Hint terminates at second order as Hi^t is a one-body 



interaction Hamiltonian. Expanding _ffi„t and dropping 
E for simplicity, the PRCC equation assumes the form 



^DCB^^W 



l^fo) 



D 

1 

'2 
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The equations of T^^-* are obtained after projecting the 
equation on singly excited states ($^|. These excitation 
states, however, must be opposite in parity to |$o)- The 
Tj^' equations are obtain in a similar way after project- 
ing on the doubly excited states ('I'^^l. The equations 
form a set of coupled nonlinear algebraic equations. The 
equations and a description of the different terms along 
with diagrammatic analysis are given in our previous 
works |11[ 112) . An approximation which incorporates all 
the important many-body effects like core-polarization is 
the linearized PRCC (LPRCC). In this approximation, 
only the terms linear in T^*^) , equivalent to retaining only 

{^/nT^i)} and {J^'^^')} in the PRCC equations. 
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FIG. 1. Diagrams of the a in the PRCC theory. The single 
excitation operators with a wavy line represent Tj^'. Simi- 
larly, the double excitation diagrams with an extra vertical 
line represent T 
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III. DIPOLE POLARIZABILITY 

In the PRCC theory we can write the a of the ground 
state of a closed-shell atom as fTTl IT2l 



a 



($o|T(i)^D + DT(i)|$o) 
(^ol^o) 



(19) 



where, D — e^'^ '^De^'^^ represents the unitary trans- 
formed electric dipole operator. Retaining only the dom- 
inant terms, we obtain 



--^($o|TPD + DTi 
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where TV = (^-ol exp[T(o)t] exp[T(o)]|$o) is the normal- 
ization factor, which involves a non-terminating series 

of contractions between T*^*^-*^ and T^^\ However, in 

the present work we use N ~ (<i>o|l + 



T2^^^T2^^\^o)- An evident advantage of computing a 
with PRCC theory is the absence of summation over the 
intermediate states |^/). The summation is subsumed 
in the evaluation of T*-^^ in a natural way and eliminates 
the need for a complete set of intermediate states. 

For further analysis and evaluation of the different 



j.^(0)ty^(0) 



terms in Eq. ( 20 ) , we use many-body diagrams or Gold- 
stone diagrams. To evaluate the diagrams we follow the 
notations and conventions given in ref. [28 . However, 
there is an additional feature in the diagrams of a, we 
employ a wavy interaction line to represent the diagrams 
of t''i \ so that it is different from the diagrams 
Similarly, to represent Tj^-* we introduce a vertical line 
to the interaction line. After due consideration of the 



equivalent diagrams, the terms in Eq. (20) correspond 
to 22 unique Goldstone diagrams and these are shown in 
Fig. [l] The equivalent algebraic expression is 
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where da 

-I-P1 

'■ab 



(a|d|6), and 
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and tZ 



^ab ^'"'^ antysymmetrized cluster amplitudes. 



In the figure, the first two diagrams. Fig. [T|a) andfljb), 
are the most important ones. These represent Tj^^D 
and DT^^\ respectively, and subsume DF and the effects 



of random phase approximation (RPA). The next two 
diagrams in the figure, Fig{I|c) and FigjIJd), arise from 
the term T^^^'^'DT^\ Similarly, the diagrams in Figllle- 

f) correspond to the hermitian conjugate, T!f'^'^'D'Ti\ 
These are the two leading order terms among the second 
order contributions, in terms of the cluster amplitudes, 
to a. The reason is, both the terms consist of dominant 
RCC and PRCC amplitudes, Tif^ and T^^\ respectively. 

Among the second order contributions, the next to 
leading order terms are T^^^'^'DTj^"^ and T!f'^^'DTf\ 
Each of these terms generate four diagrams, which are 
given in Fig{lJo-r) and Figjljs-v) and these correspond to 
T^^^'^DTj^"^ and Tlf'^^'DT^^, respectively. The remaining 
second order terms, T^^'^'^ 'Dt[°\ T^^^'^'DtI^^ and their 
hermitian conjugates, have marginal contributions to a. 
However, for completeness, these are included in the com- 
putations. 



TABLE I. The ao and P parameters of the even tempered 
GTO basis used in the present calculations. 



Atom 


s 




P 




d 










QO 


P 


ceo 




Na+ 


0.0025 


2.210 


0.00955 


2.125 


0.00700 


2.750 


K+ 


0.0055 


2.250 


0.00995 


2.155 


0.00690 


2.550 


Rb+ 


0.0052 


2.300 


0.00855 


2.205 


0.00750 


2.145 


Cs+ 


0.0097 


2.050 


0.00975 


2.005 


0.00995 


1.705 


Fr+ 


0.0068 


2.110 


0.00645 


2.050 


0.00985 


1.915 



IV. CALCULATIONAL DETAILS 

A. Basis set and nuclear density 

The first step of our computations, which is also true 
of any atomic and molecular computations, is to gener- 
ate an spin-orbital basis set. For the present work, the 
basis set is even-tempered gaussian type orbitals (GTOs) 
P5] generated with the Dirac-Hartree-Fock Hamiltonian. 
This means, the radial part of the spin-orbitals are linear 
combinations of the Gaussian type functions. The Gaus- 
sian type functions which constitutes the large compo- 
nents are of the form 



^Kp 



(22) 



where p = 0, 1 ... m is the GTO index and m is the 
number of gaussian type functions. The exponent ap — 
aol3P~^, where ao and /3 are two independent paramters. 
Similarly, the small components of the spin-orbitals are 
linear combination of g^p{r), which are generated from 
g^p(r) through the kinetic balance condition [5D]. The 
GTOs are calculated on a grid [3T] with optimized val- 
ues of Qfo and /3. The optimization is done for individual 
atoms to match the spin-orbital energies and self con- 
sistent field (SCF) energy of GRASP92 [32 . For the 
current work, the optimized ao and /3 are listed in Table. 
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Hj For comparison, the spin-orbital energies of Cs"*" ob- 
tained from the GTO and GRASP92 are hsted in Table 
|TTj In the table, the deviation of the GTO results from 
the GRASP92 is 10"^, which is quite smah. We obtain 
similar level of deviations for the other ions as well. 



TABLE II. Core orbital energies of Cs^ in atomic units. 



Orbital 


DC 


GRASP92 32 


lSl/2 


-1330.1173 


-1330.1129 


2Si/2 


-212.5643 


-212.5673 


2pi/2 


-199.4294 


-199.4288 


2p3/2 


-186.4366 


-186.4358 


3Si/2 


-45.9697 


-45.9695 


3pi/2 


-40.4483 


-40.4455 


3p3/2 


-37.8943 


-37.8917 


3^3/2 


-28.3096 


-28.3030 


3^5/2 


-27.7752 


-27.7682 


4Si/2 


-9.5128 


-9.5106 


4pi/2 


-7.4463 


-7.4437 


4p3/2 


-6.9209 


-6.9188 


4^3/2 


-3.4856 


-3.4921 


4^5/2 


-3.3969 


-3.4038 


5Si/2 


-1.4898 


-1.4933 


5pi/2 


-0.9079 


-0.9139 


5p3/2 


-0.8403 


-0.8459 



The next step, related to spin-orbial basis set, is the 
choice of an ideal basis set size. For this, we examine the 
convergence of a using the LPRCC theory. We calculate 
a starting with a basis set of 50 GTOs and increase the 
basis set size in steps through a series of calculations. 
The results of the such a series of calculations are listed 



in Table. Ill it shows the convergence of a of Cs"*" as a 
function of basis set size. 

In the present work we have considered finite size Fermi 
density distribution of the nucleus 



Po 



1 -I- e(''-c)/a ' 



(23) 



where, a = <41n(3). The parameter c is the half charge 
radius so that Pnuc(c) — Pq/2 and t is the skin thickness. 
Coming to the PRCC equations, these are solved itera- 
tively using Jacobi method, we have chosen this method 
as it is parallelizable. The method, however, is slow to 
converge. So, to accelerate the convergence we use direct 
inversion in the iterated subspace (DIIS) [SUI- 



TABLE III. Convergence pattern of a (Cs"*") as a function of 
the basis set size. 



No. of orbitals 



Basis size 



103 
117 
131 
143 
159 
177 



(15s, 13p, 13d,9/,9g) 14.9480 

(17s, 15p, 15d, 11 f, 9g) 14.9235 

(19s,17p, 17d,ll/,ll5) 14.9124 

(23s,19p, 19d,13/,llff) 14.9086 

(23s,21p,21d,13/,13ff) 14.9077 

(25s, 23p, 23d, 15/, 15g) 14.9077 



B. Intermediate Diagrams 

The PRCC diagrams corresponding to the nonlinear 
terms are numerous and topologically complex. Fur- 
ther more, in these diagrams, the number of the spin- 
orbitals involved is large and in general, the diagrams 
with the largest number of spin-orbitals are associ- 



ated with the terms HnT!;°''t!-°'> 



_^^y^(0) j,^(0) y|0) rp{0) 



All of these terms have a com- 
mon feature: the presence of the Coulomb integral 
{ab\l/ri2\pq) ■ Returning to the number of spin-orbitals, 
the diagrams arising from any of the three terms 
mentioned earlier consist of four core and virtual spin- 
orbitals each. Accordingly, the number of times a dia- 
gram is evaluated, N^, scales as n'^nf, and sets the com- 
putational requirements. Here, Uo and n„ are the num- 
ber of core and virtual spin-orbitals, respectively. In the 
present work, Uq ^ 10 and n„ ^ 100 for lighter atoms and 
moderate sized basis sets, even then ~ 10^^. This is 
a large number and puts a huge constraint on the com- 
putational resources. 





FIG. 2. 



Example diagrams of Hj^T^^T^^ 



which contribute 

to the Tj^' equations. The portion of the diagrams within 
the rectangles with rounded corners are examples of the one 
body (a) and two-body (b) intermediate diagrams. 



To mitigate the computational constraints arising from 
the n^n^ scaling, we separate the diagrams into two 
parts. One of the parts scales at the most n^n^ and the 
total diagram is equivalent to the product of the parts. 
The part of the diagram, which is calculated first is re- 
ferred to as the intermediate diagram. During computa- 
tions, all the intermediate diagrams are calculated first 
and stored. Later, these are used to combine with the re- 
maining part of the RCC diagram and the total diagram 
is calculated. The scaling is still n^n^ and compared to 
the n^n^ scalling, this improves the performance by sev- 
eral orders of magnitudes. 

To examine in more detail, consider the term 
i?jv72°^T2^-', the algebraic expression for one of the terms 

.(0) 



contributing to the is 



i'^ab)2e,"'i>"'a'^baa 



pq t 



(24) 



and it is diagrammatically equivalent to Fig. [2|^a). How- 
ever, while evaluating the diagram, the part within the 
dashed round rectangle or the intermediate diagram can 
be separated and computed first. The Eq. (24) can then 
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be written as 



2a O'Wq'^baa 



E 

d 



(25) 



where ryf = J2rcs ^acKs is the amphtude of the effective 
one-body operator corresponding to the intermediate di- 
agram. It scales as n^n^ and when contracted with 
the computation stih cales as n^n^- This is much less 
than the n'^n'^ scaling. Consider another term 



("^ab )2b (^Wq(^f>aa 



T.^alvfycH4^baa, (26) 



and it is diagrammatically equivalent to Fig. ^h). Like 
in the previous case, the intermediate diagram ( part 
within the dashed round-rectangle ) can be calculated 
first and the equation can rewritten as 



cd 



(27) 

Here, the intermediate diagram corresponds to a two- 



body effective operator with amphtude ?7^^ = J2rs ^al^rs 
and scales as n^n^- The scaling remains the same when 
the total diagram is evaluated. Extending the method 
to other diagrams, depending on the topology, there are 
other forms of one-body and two-body intermediate dia- 
grams. 




FIG. 3. Contributions to the next to leading order terms 
Ti'^'^DTj"' + h.c. in terms of the pairs of core spin-orbitals. 



V. RESULTS AND DISCUSSIONS 

To compute a using PRCC theory, as described earlier, 
we consider terms up to second order in the cluster oper- 
ators. We have, however, studied terms which are third 
order in cluster operators and examined the contribu- 
tions from the leading order terms. But the contributions 



are negligible and this validates our choice of considering 
terms only upto second order in cluster operators. To 
begin with, we compute a using the cluster amplitude 
obtained from the LPRCC and results are presented in 



Table IV In the table we have listed, for systematic com- 
parison, the experimental data and results from previous 
theoretical computations. 

TABLE IV. Static dipole polarizabilty of alkali ions 



Atom LPRCC RCCSDT 34 RRPA 35 



Expt. 



Na+ 


1.009 


1.00(4) 


0.9457 


0.9980(33) 


K+ 


5.521 


5.52(4) 


5.457 


5.47(5) [37] 


Rb+ 


8.986 


9.11(4) 


9.076 


9.0 [38] 


Cs+ 


14.924 


15.8(1) 


15.81 


15.644(5) [39] 


Fr+ 


19.506 


20.4(2) 







For Na^ and K+, our values of a are higher than the 
experimental values by 1% and 0.9% respectively. How- 
ever, for Rb+ and Cs+ our results are lower than the 
experimental values by 0.15% and 4.8%, respectively. In 
terms of theoretical results, our results of Na"*" and K+ 
are in excellent agreement with the previous work which 
used the RCCSDT method for computation. But, for 
Rb'^ and Cs^, like in the experimental data, our results 
are lower than the RCCSDT results. One possible reason 
for these deviations in the heavier ions could be the exclu- 
sion of triple excitation cluster operators in the present 
work. Our result of Fr+ seems to bear out this reasoning 
as the same trend is observed ( our result is 4.4% lower 
than the RCSSDT result) in this case as well. However, 
in absence of experimental data for Fr+, it is difficult to 
arrive at a definite conclusion. 

To investigate the importance of Breit interaction, , 
in computing a of the alkali ions, we exclude in the 
atomic Hamiltonian and do a set of systematic calcula- 
tions. Our results for the values of a are then 1.008, 
5.514, 8.973 and 14.908 for Na+, K+, Rb+ and Cs+, re- 
spectively These values are 0.001, 0.007, 0.013 and 0.016 
a.u lower than the results computed using the Dirac- 
Coulomb-Breit Hamiltonian. This indicates that the cor- 
rection from the Breit interaction is larger in heavier ions 
and this is as expected since the stronger nuclear poten- 
tial in heavier ions translates to larger Breit correction. 

For a more detailed study, we examine the contribu- 



tions from each of the terms in the Eq. 20 and these are 
listed in Table. |V] The leading order contribution arises 
from T^^-'^D -|- h.c and diagrammaticaly, it corresponds 
to the first two diagrams in Fig. [1] . These are also 
the lowest order terms and are the dominant terms since 
these subsume the contributions from the Dirac-Fock and 
RPA effects. For all the ions, the results from the dom- 
inant terms exceeds the final results. Here, it must be 
mentioned that a similar trend is observed in the results 
of noble gas atoms as well [TT] [T^] . The next to leading 

order (NLO) contributions arise from the Ti^^)^TiT!f^ 
+ h.c. The contributions from these terms are an or- 
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TABLE V. Contribution to a from different terms and their 
hermitian conjugates in the hnearized PRCC theory. 



.(i)t 



D + h.c. and The dominant contributions are pre- 
sented in Table. IVIII It is to be noted that a has 



Terms + h.c. 



"Rb^ 



Cs^ 



Fr^ 





1.018 


6.043 


10.029 


17.472 


22.926 


T^(l)tj3yTO 


-0.018 


-0.299 


-0.519 


-1.023 


-1.326 




0.012 


-0.038 


-0.072 


-0.188 


-0.126 




-0.001 


0.008 


0.016 


0.039 


0.026 




0.023 


0.204 


0.332 


0.654 


0.834 


Normahzation 


1.025 


1.072 


1.089 


1.136 


1.145 


Total 


1.009 


5.521 


8.986 


14.924 


19.506 



der of magnitude smaller then T^^'^^^D + h.c but more 
importantly, the contributions are opposite in phase. In- 
terestingly, the next important terms T2 (i)tDr(°) + h.c 
have contributions which nearly cancels the NLO con- 
tributions. Continuing further, among the second order 
terms, the smallest contribution arise from Ts^i^^DTi'"^ 
+ h.c, which is perhaps not surprising since T2^^-'^ and 
T^^^ are the cluster operators with smaller amplitudes in 
PRCC and RCC theories, respectively. Collecting the re- 
sults, the net contributions from the second order terms 
are 0.016, -0.117, -0.223, -0.456 and -0.517 for Na+, K+, 
Rb^, Cs^ and Fr"*", respectively. 

Next, we consider all the terms in the PRCC theory, 
including the terms which are non-linear in cluster op- 
erators. The results of a are presented in the Table |VI| 
For Na^ the result of a is 2.6% higher than the exper- 



TABLE VI. Contribution to a from different terms and their 
conjugate in the PRCC theory after including the terms non- 
linear in cluster operators. 



Terms + h.c. 



Na-t 



Ic^ rF^ 



"Ci^ 





1.034 


6.302 


10.438 


18.376 




-0.018 


-0.316 


-0.544 


-1.084 




0.012 


-0.040 


-0.076 


-0.198 




-0.0008 


0.008 


0.016 


0.038 




0.024 


0.194 


0.308 


0.596 


Normalization 


1.026 


1.072 


1.090 


1.136 


Total 


1.025 


5.735 


9.305 


15.606 



imental value. Similarly, for K"*" and Rb^ the results 
are 4.6% and 3.3% higher than the experimental values. 
For Cs^ the nonlinear PRCC theory gives a much im- 
proved result than the LPRCC results and the deviation 
from the experimental value is reduced to 0.24%. On a 
closer examination, most of the change associated with 
the nonlinear PRCC can be attributed to the increased 
contribution from T^^-'^D -I- h.c. As these terms sub- 
sume RPA effects, the increased contributions indicate 
that RPA effects are larger in the nonlinear PRCC. 

To investigate the RPA effects in detail, we isolate 
the contributions from each of the core spin-orbitals to 



TABLE VII. Four leading contributions to {T^^'^D + h.c} to 
a in terms of the core spin-orbitals. 



Na+ 


K+ 


Rb+ 


0.652 (2p3/2) 
0.322 (2pi/2) 
0.044 (2S1/2) 
0.0004 (lsi/2) 


4.016 (3p3/2) 
1.938 (3pi/2) 
0.076 (3si/2) 
0.008 (2p3/2) 


6.858 (4p3/2) 
3.038 (4pi/2) 
0.058 (4si/2) 
0.044 (3d5/2) 


Cs+ 


Fr+ 




12.375 (5p3/2) 
4.735 (5pi/2) 
0.192 (4^5/2) 
0.121 (4^3/2) 


18.287 (6P3/2) 
4.073 (6P1/2) 
0.376 (5^5/2) 
0.211 (5d3/2) 





a quadratic dependence on the radial distance, so the 
orbitals with larger spatial extension contribute domi- 
nantly. The effect of this is discernible in the results, for 
all the alkali ions the leading contribution to a comes 
from the outermost orbital, which is the occupied 

orbital with largest radial extent. The next leading con- 
tribution arise from the npi/2 orbital. An important ob- 
servation is, as we proceed from from lower Z to higher Z, 
the ratio of contribution of 71.^3/2 to the npi/2 increases. 
It is 1.8, 2.1, 2.3, 2.6 and 4.5 for Na+, K+, Rb+, Cs+ 
and Fr^ respectively. The ratio is much larger in the 
case Fr^ and without any ambiguity it can be attributed 
to the relativistic contraction of the npi/2 orbital. The 
third leading contribution for Na^, K+, Rb^ arise from 
the 2si/2, 3si/2 and 4sx/2 orbital respectively. But, for 
Cs"*" and Fr"*" the third leading contribution arise from 
4(^5/2 and 5^5/2 orbital respectively. This is because the 
5si/2 and 6si/2 orbital are contracted due to large rel- 
ativistic effects. From the above analysis of RPA effets, 
the trend in the contributions demonstrates the impor- 
tance of relativistic corrections in Cs"*" and Fr"*". 



TABLE VIII. Core orbitals contribution from Ti(^^t£)j.^(0) 
a of Na+ and K+ 



Na+ 


K+ 


-0.0040 (2p3/2,2p3/2) 
-0.0021 (2p3/2,2pi/2) 
-0.0021 (2pi/2,2p3/2) 
-0.0010 (2pi/2,2pi/2) 


-0.0646 (3p3/2,3p3/2) 
-0.0367 (3p3/2,3pi/2) 
-0.0360 (3pi/2,3p3/2) 
-0.0148 (3pi/2,3pi/2) 



To study the pair-correlation effects, we identify the 
pairs of core spin-orbitals in the next leading order terms 
Ti(i)1'DT2^°) + h.c. The four leading order p airs fo r Na 
and K+, Rb+ 



VIII 



and 



IX 



Cs^ and Fr^ are listed in table 
respectively. The dominant contribution, for all the ions, 
arise from the combination (np3/2, '^^3/2) orbital pairing. 
To illustrate the relative values, the contributions from 



8 



the pairs of the five outermost core spin-orbitals of Rb+ is 



shown as a barchart in Fig. IV B Comparing the results 



TABLE IX. Core orbitals contribution from Ti'-'^^^DT^"' to 
a of Rb+, Cs+ and Fr+ 



Rb+ 


Cs+ 


-0.1113 (4p3/2,4p3/2) 

-0.0601 (4p3/2,4pi/2) 
-0.0565 (4pi/2,4p3/2) 
-0.0223 (4pi/2,4pi/2) 


-0.2126 (5p3/2,5p3/2) 
-0.1073 (5p3/2,5pi/2) 
-0.0930 (5pi/2,5p3/2) 
-0.0347 (5pi/2,5pi/2) 


Fr+ 


-0.3078 (6p3/2,6p3/2) 
-0.1266 (6p3/2,6pi/2) 
-0.0828 (6pi/2,6p3/2) 
-0.0489 (6p3/2,5d5/2) 



of all the ions, there is a major difference in the results 
of Fr+. For Fr+ the fourth largest contribution is from 
the (6p3/2, 5^5/2) pair, whereas for the other ions it is of 
the form 't-Pi/2)- This is again a consequence of 

the contraction of the 6si/2 spin-orbital in Fr+ due to 
relativistic effects. 

In the present calculations we have identified the fol- 
lowing possible sources of uncertainty. The truncation of 
the spin-orbital basis sets is one of the possible source. 
For all the ions we start the computations with 9 symme- 
tries and increase up to 13 symmetries. Along with it, we 
also vary the number of the spin-orbitals till a converges 
to sa 10^''. So, we can safely neglect this uncertainty 
for our calculations. Another source of uncertainty is 
the truncation of the CC theory at the single and dou- 
ble excitation for both unperturbed and the perturbed 
RCC theories. Based on our previous theoretical results 
PH, [T^ the contributions from triple and higher order 
excitations is at the most ~ 3.3%. The truncation of 
De +e De 



' Ti^^* ' I ^t'"' ti^t' ' -j-j^g second order in cluster 



operator is also a source of uncertainty. From our ear- 
lier studies [26] with CC theory and in the present work 
we have studied the contribution from the third order in 
cluster operator, but the contribution is negligibly small. 
The quantum electrodynamical(QED) corrections is an- 
other source of uncertainty in our calculations and based 
on our previous studies, we estimate it at 0.1%. In total, 
we estimate the uncertainties in our results as ~3.4%. 



VI. CONCLUSION 

We have computed the static electric dipole polariz- 
ability of alkali ions using the PRCC theory. The PRCC 
theory is a coupled-cluster based theory and can be eas- 
ily modified to incorporate other perturbations in the 
atomic many-body calculations. In the present work, we 
have explored the use of PRCC theory to calculate the 
electric dipole polarizability of closed-shell ions and find 
that the results are in good agreement with the experi- 
mental results and previous theoretical results. 

On a closer examination of the results, the pattern 
of the contributions from the individual and pairs of 
spin-orbitals establishes the importance of the relativis- 
tic corrections in higher Z ions. The results further in- 
dicates that it is essential to obtain the outermost P3/2 
spin-orbitals of the ions accurately. The reason is, these 
are associated with the dominant contributions from the 
Dirac-Fock, RPA effects and pair-correlation effects. 
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